The correlation structure of spatial autoregressions

This paper investigates how the correlations implied by a first-order simultaneous autoregressive (SAR(1)) process are affected by the weights matrix and the autocorrelation parameter. A graph theoretic representation of the covariances in terms of walks connecting the spatial units helps to clarify a number of correlation properties of the processes. In particular, we study some implications of row-standardizing the weights matrix, the dependence of the correlations on graph distance, and the behavior of the correlations at the extremes of the parameter space. Throughout the analysis differences between directed and undirected networks are emphasized. The graph theoretic representation also clarifies why it is difficult to relate properties ofW to correlation properties of SAR(1) models defined on irregular lattices

Testing for spatial autocorrelation: the regressors that make the power disappear

We show that for any sample size, any size of the test, and any weights matrix outside a small class of exceptions, there exists a positive measure set of regression spaces such that the power of the Cliff-Ord test vanishes as the autocorrelation increases in a spatial error model. This result extends to the tests that define the Gaussian power envelope of all invariant tests for residual spatial autocorrelation. In most cases, the regression spaces such that the problem occurs depend on the size of the test, but there also exist regression spaces such that the power vanishes regardless of the size. A characterization of such particularly hostile regression spaces is provided.Cliff-Ord test; point optimal tests; power; spatial error model; spatial lag model; spatial unit root

Spatial design matrices and associated quadratic forms: structure and properties

The paper provides significant simplifications and extensions of results obtained by Gorsich, Genton, and Strang (J. Multivariate Anal. 80 (2002) 138) on the structure of spatial design matrices. These are the matrices implicitly defined by quadratic forms that arise naturally in modelling intrinsically stationary and isotropic spatial processes. We give concise structural formulae for these matrices, and simple generating functions for them. The generating functions provide formulae for the cumulants of the quadratic forms of interest when the process is Gaussian, second-order stationary and isotropic. We use these to study the statistical properties of the associated quadratic forms, in particular those of the classical variogram estimator, under several assumptions about the actual variogram.Cumulant, Intrinsically Stationary Process, Kronecker

Spatial design matrices and associated quadratic forms: structure and properties

The paper provides significant simplifications and extensions of results obtained by Gorsich, Genton, and Strang (J. Multivariate Anal. 80 (2002) 138) on the structure of spatial design matrices. These are the matrices implicitly defined by quadratic forms that arise naturally in modelling intrinsically stationary and isotropic spatial processes.We give concise structural formulae for these matrices, and simple generating functions for them. The generating functions provide formulae for the cumulants of the quadratic forms of interest when the process is Gaussian, second-order stationary and isotropic. We use these to study the statistical properties of the associated quadratic forms, in particular those of the classical variogram estimator, under several assumptions about the actual variogram.Cumulant; Intrinsically stationary process; Kronecker product; Quadratic form; Spatial design matrix; Variogram

Spatial circular matrices, with applications

The cumulants of the quadratic forms associated to the so-called spatial design matrices are often needed for inference in the context of isotropic processes on uniform grids. Unfortunately, because the eigenvalues of the matrices involved are generally unknown, the computation of the cumulants may be very demanding if the grids are large. This paper constructs circular counterparts, with known eigenvalues, to the spatial design matrices. It then studies some of their properties, and analyzes their performance in a number of applications.

Non-Identifiability in Network Autoregressions

We study identification in autoregressions defined on a general network. Most identification conditions that are available for these models either rely on repeated observations, are only sufficient, or require strong distributional assumptions. We derive conditions that apply even if only one observation of a network is available, are necessary and sufficient for identification, and require weak distributional assumptions. We find that the models are generically identified even without repeated observations, and analyze the combinations of the interaction matrix and the regressor matrix for which identification fails. This is done both in the original model and after certain transformations in the sample space, the latter case being important for some fixed effects specifications

Non-testability of equal weights spatial dependence

We show that any invariant test for spatial autocorrelation in a spatial error or spatial lag model with equal weights matrix has power equal to size. This result holds under the assumption of an elliptical distribution. Under Gaussianity, we also show that any test whose power is larger than its size for at least one point in the parameter space must be biased

Exact and higher-order properties of the MLE in spatial autoregressive models, with applications to inference.

The quasi-maximum likelihood estimator for the autoregressive parameter in a spatial autoregression usually cannot be written explicitly in terms of the data. A rigorous analysis of the first-order asymptotic properties of the estimator, under some assumptions on the evolution of the spatial design matrix, is available in Lee (2004), but very little is known about its exact or higher-order properties. In this paper we first show that the exact cumulative distribution function of the estimator can, under mild assumptions, be written in terms of that of a particular quadratic form. Simple examples are used to illustrate important exact properties of the estimator that follow from this representation. In general models a complete exact analysis is not possible, but a higher-order (saddlepoint) approximation is made available by the main result. We use this approximation to construct confidence intervals for the autoregressive parameter. Coverage properties of the proposed confidence intervals are studied by Monte Carlo simulation, and are found to be excellent in a variety of circumstance